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We studied the electronic structure Bi 2 Se 3 employing density functional theory. The calculations 
show that the Dirac states primarily consists of the states at the interface of surface and sub-surface 
quintuple layers and the emergence of the Dirac states depends on the surface terminations in sharp 
contrast to their surface character expected for such systems. This manifests the complexity of real 
materials and the realization of topological order in real materials as an outstanding problem. We 
discover that the surface character of the Dirac states can be achieved by adsorption of oxygen on 
Bi-terminated surface due to the change in covalency by the relatively more electronegative oxygens. 

PACS numbers: 73.20.At, 03.65.Vf, 68.43.-h 


Topological insulators are characterized by metallic 
states at the surface of bulk insulators protected by time 
reversal symmetry and strong spin-orbit coupling [Ei- 
Due to time reversal symmetry protection, the electrons 
moving on the surface cannot be scattered in the re¬ 
verse direction without spin flip. Thus, the corresponding 
energy bands exhibit linear dispersion forming a Dirac 
cone, which is one of the key signatures in the electronic 
structure of a topological insulator. Due to such unique 
electronic properties, these materials are potential candi¬ 
dates for realization of exotic physics such as Majorana 
Fermions, magnetic monopoles [31,1^. These materials are 
also expected to bring immense technological advances 
and new possibilities in the fields of spintronics, quan¬ 
tum computation, dissipation less charge transfer etc. [H . 
While all these realizations require experimental finding 
of protected topologically ordered surface states on insu¬ 
lating bulk materials, experiments show contrasting sce¬ 
narios often with instability of topological order. The 
bulk of almost all the materials studied exhibit metallic 
phase [H. Thus, finding of a true topological insulator 
remains to be an outstanding puzzle. 


(a) Crystal structure 


(b) Surface and bulk BZ 



FIG. 1. (a) Crystal structure of Bi 2 Se 3 exhibiting atomic ar¬ 
rangement and the rhombohedral unit cell. SI, S2, S3 and 
S4 denote the layers where the surface terminations are con¬ 
sidered in this work, (b) The surface and bulk Brillouin zone 
(BZ) along with different high symmetry points. 


Engineering new materials for such behavior requires 
critical understanding of the electronic structure of the 
materials exhibiting topological order. For example, one 
of the most studied materials for such behavior is Bi 2 Se 3 
that shows Dirac cone with its apex at finite binding en¬ 
ergy instead of its chemical potential. In addition, sev¬ 
eral bulk bands are observed at the Fermi level suggest¬ 
ing its metallic bulk electronic structure. Carrier doping 
due to aging/impurities often inhibits achieving topologi¬ 
cal transport regime. Some theories predicted relaxation 
of van der Waals bond as a reason for such behavior 
while other studies show necessity of unusually large 
change in bond length to produce the experimentally ob¬ 
served scenario Q . The surface band bending induced by 
adsorption of residual gases , commonly observed in 
semiconductors, may also give rise to such effect. Inter¬ 
estingly, both positive and negative charge doping due to 
impurity adsorption have been observed in angle resolved 


photoemission measurements 11 fH- While all these 
studies indicate instability of the surface states with ag¬ 
ing, impurity adsorption etc., there are studies reporting 
stable surface behavior |14l4l6|. Evidently, the behavior 
of the surface states of these materials is highly anoma¬ 
lous. In order to reveal microscopic details underlying 
such electronic properties, we calculated the electronic 
band structure for both surface and bulk employing den¬ 
sity functional theory, and discover interesting results, 
significantly different from the expected behavior of this 
system. 


The electronic band structure calculation of Bi 2 Se 3 
was carried out using full potential linearized augmented 
plane wave method within the local density approxima¬ 
tions EE- The bulk electronic structure is calculated 
using the lattice constants, a = 5 = 4.18 A, c = 28.7 A, 
a = l3 = 90*^ and 7 = 120 '^ obtained from the literature 
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FIG. 2. Band structure of (a) clean and (b) 1 monolayer (ML) O covered SI terminated surfaces, and (c) clean and (d) 1 ML 
O covered S2 terminated surface. The black dots show the bands with no spin-orbit coupling and lines are with spin-orbit 
coupling. The blue region is the projected bulk bands and the red dots denote the bands with more than 50% O contribution. 


[rij . In order to investigate the surface electronic struc¬ 
ture, we considered the lattice in a ’slab’ configuration 
with at least 5 quintuple layers within the unit cell - the 
total number of layers varies for different surface termi¬ 
nations. Calculations were carried out for all the possible 
surface terminations as defined in Fig. 1. 

Bi 2 Se 3 forms in rhombohedral structure with space 
group, R3m. In Fig. 1(a), we show the crystal 
structure [lil with both hexagonal and rhombohedral 
axes. It has a layered structure containing quintuple 
layers consisting of Sel-Bi-Se2-Bi-Sel layers, where two 
nonequivalent Se layers are denoted by Sel and Se2. The 
quintuple layers are believed to be connected through 
weak Van der Waals bond the inter-layer bond¬ 

ing within the quintuple layer is relatively stronger. This 
indicates a well defined cleavage plane exposing Sel ter¬ 
minated surface. Recent experiments, however, show sig- 
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The calculations for the bulk Bi 2 Se 3 converge to in¬ 
sulating ground state with a band gap of about 0.3 eV 
consistent with the experimental findings 23-25|. The 
results from the slab calculations exhibit interesting sce¬ 
nario with very different surface bands near the Fermi 
level, ep for different terminations. All the energy bands 
obtained from slab calculations are doubly degenerate as 
a symmetric slab geometry has been used [26| for these 
calculations. In Figs. 2(a) and 2(c), we show the band 
structure obtained for SI and S2 cases. The black dots 
represent the bands without spin-orbit coupling (SOC) 
and the solid red lines correspond to the calculations in¬ 


nature of both Bi and Se terminated surfaces on cleaving 
[ni0 exhibiting significantly different behavior. Thus, 
we considered the surface terminations SI, S2, S3 and 
S4 (see figure) representing Se2, Sel (one Sel on top of 
the quintuple layer), Sel (normal cleavage plane) and Bi 
terminations, respectively. Fig. 1(b) shows the bulk and 
surface Brillouin zones with different high symmetry k 
points. For (111) surface of Bi 2 Se 3 , the time reversal in¬ 
variant momenta (TRIM) are the high symmetry points, 
F and M. In addition to the pristine cases, we deposited 
a monolayer of oxygen on top of the surfaces. The ex¬ 
change correlation potential is included under generalized 
gradient approximation (GGA) 0. Spin-orbit interac¬ 
tion is included as a second order perturbation for Bi 
and Se atoms [22| . The energy convergence was achieved 
using 10 X 10 X 1 /c-mesh. 


r 


eluding SOC. The shaded area represents the projected 
bulk bands (PBB) for the SOC case, where the difference 
in Fermi level between the surface and bulk calculations 
are compensated by matching the bulk bands away from 
the Fermi level. The bands appearing in the gap of the 
PBB are the surface bands. It is evident from the data 
that SOC leads to a splitting of the spin degenerate bands 
obtained without SOC [22|. For both the surface termi¬ 
nations, the surface bands cross the Fermi level between 
surface projected TRIMs, F-M odd number of times im¬ 
plying that the material belong to a strong topological 
insulator class Q. The bands near ep consist primarily 
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FIG. 3. Energy band structure for (a) S3-terminated, (b) 1 ML oxygen on S3 surface, (c) S4 terminated surface and (d) 1 ML 
oxygen on S4 surface, (e) Shift of slab Fermi level w.r.t. bulk Fermi level (dotted line) binding energy. 


of Se Ap and Bi 6p characters. 

For clean Se2 terminated surface, the conduction band 
maxima and the valence band minima at F point shown 
by arrows in Fig. 2(a) exhibit a separation of about 0.3 
eV due to Bi p and Se2 p hybridizations. The band cross¬ 
ing at around -0.16 eV exhibits a scenario akin to a Dirac 
cone representing topological order. Here, the finite en¬ 
ergy of the Dirac point (apex of the cone) is intrinsic to 
the electronic structure. Impurities/defects play impor¬ 
tant role in the electronic structure leading to charge car¬ 
rier doping and introducing local character of the charge 
carrier [271, l28| . We have verified this in the present case 
by calculating electronic structure with one monolayer 
oxygen on the surface. The Dirac point shifts away from 
Ci? as shown by an arrow in Fig. 2(b). An additional 
pair of bands appear near ep possessing primarily O 2p 
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The energy bands corresponding to S3 and S4 termi¬ 
nated surfaces are shown in Fig. 3. The results obtained 
without SOC exhibit finite gap between the valence and 
conduction bands. Inclusion of SOC leads to the forma¬ 
tion of distinct Dirac cone like feature at F point at 
as shown by arrow in the figure. In these terminations, 
the presence of spin-orbit coupling leads to the surface 
bands crossing the Fermi level odd number of times be¬ 
tween TRIMs, F and M providing again the signature 
of a strong topological insulating phase in these cases 
US liS. The Dirac point (DP) appears very close to 
the top of the valence band for Se terminated surface 
(S3 case) and near bottom of the conduction band on 
Bi-terminated surface (S4 case). 


character. This is shown by big circles representing more 
than 50% O 2p contributions in Fig. 2(b). The energy 
bands near ep becomes narrower compared to the pris¬ 
tine case suggesting enhancement of local character of the 
electronic states. The overall shift of the energy bands 
towards higher binding energies (= negative of the energy 
shown in the figures) in the later case represents an effec¬ 
tive electron doping as expected due to relatively higher 
electronegativity of oxygen. The results for S2 termina¬ 
tion are shown in Fig. 2(c) and 2(d). The energy bands 
exhibit significantly different behavior with no signature 
of Dirac cone. Oxygen in this case leads to a band nar¬ 
rowing (see Fig. 2(d)) and effective electron doping as 
found in the case of SI termination. However, the oxy¬ 
gen 2p character appears stronger around high symmetry 
points other than F. 


The most interesting phenomena occurs on deposition 
of oxygen on these surfaces; the Dirac cone shifts to 
higher binding energy (^ 0.14 eV) along with signifi¬ 
cant narrowing of the energy bands near ep (see Fig. 
3(b)). O 2p-derived additional narrow bands appear in 
the vicinity of the Fermi level as found in the case of Se2 
terminated surface. Fig. 3(d) exhibits the energy bands 
corresponding to a monolayer oxygen on Bi-terminated 
surface. The DP shifts away from the bottom of the con¬ 
duction band and the Dirac cone becomes prominent, 
better defined with larger slope. Here, the oxygen 2p- 
derived bands with smaller dispersion appear below ep 
unlike the Se-terminated cases. This makes the topolog¬ 
ical surface states better accessible and easily discernible 
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FIG. 4. Angle resolved photoemission (ARPES) data adapted 
from ref. [ill] . ARPES data from Se-terminated surface (a) 
freshly cleaved and (b) 40 hours aged, and Bi-terminated sur¬ 
face (c) freshly cleaved and (d) 40 hours aged. Inset shows 
Dirac cone region with magnified intensity. 


from the impurity bands. 

In Fig. 3(e), we show the Fermi energy found in the 
slab calculations relative to that in the bulk calculation. 
While SI, S2 and S4 cases show an increase in Fermi 
energy, the Fermi energy becomes smaller than the bulk 
Fermi energy in S3 case. The S4 case (Bi terminated 
surface) exhibits the highest energy difference. Presence 
of 1 ML oxygen leads to a shift of the Fermi level towards 
higher energies suggesting an effective electron doping 
into the system for all the Se terminations. Ironically, the 
deposition of oxygen on Bi terminated surface (S4 case) 
leads to a shift of the Fermi level in opposite direction - 
towards lower energies suggesting an effective hole doping 
in this case consistent with the experimental observations 
0 . 

It is clear that the electronic structure corresponding 
to different surface terminations of the system exhibit 
different scenario. The sample with Sel bi-layer on top 
(S2 case) do not show Dirac cone although there are odd 
number of fermi level crossing of the energy bands be¬ 
tween TRIMs. Here the top Se layers are weakly coupled 
with the bulk quintuple layer underneath and the surface 
electronic structure possess effectively two dimensional 
nature. 

The samples with Se2, Sel and Bi terminated surface 
as discussed in SI, S3 and S4 cases exhibit Dirac cone 
in their energy band structure. The presence of a layer 
of oxygen leads to instability of the Dirac states corre¬ 
sponding to the Se-terminated surface. In all the cases, 
the surface oxygen bands appear close to the Fermi level 
along with a shift of the Dirac cone to higher binding 
energies 0[i3. Interestingly, the Bi-terminated surface 
exhibit more prominent Dirac states on oxygen deposited 
surface. In order to investigate the correspondence of this 
finding with the experimental scenario, we show the ex¬ 
perimental angle resolved photoemission results in Fig. 
4 adapted from the Ref.[ll|. The experimental results 



FIG. 5. Gontribution to the Dirac point from different atomic 
layers, (a) S3, (b) S4 and (c) O on S4 terminations. 


for the Se-terminated (Clvl) and Bi-terminated (Clv2) 
cases are shown in upper and lower panels of Fig. 4, re¬ 
spectively. The Freshly cleaved surface exhibit signature 
of Dirac cone in both the cases as expected. Oxygen ac¬ 
cumulation on the Se-terminated surface leads to a shift 
of the Dirac point towards higher binding energy due 
to effective electron doping as found in this study too, 
along with a weakly defined Dirac cone. However, the 
Bi-terminated surface exhibits energy shift in opposite 
direction as predicted in Fig. 3(e) and exhibits a better 
defined conical shape of the Dirac cone. These results 
clearly corroborate the theoretical description shown in 
Fig. 3 establishing the robust nature of the Dirac states 
on Bi-terminated surface. 

In order to understand the above scenario better, we 
investigate the character of the Dirac states in Fig. 5. 
The conduction band in this system consists of Se 4p and 
Bi 6p states. The Se 4p and Bi 6p contributions to the 
Dirac point (intensity within ±50 meV across the Dirac 
point) obtained for different layers of Sel terminated sur¬ 
face is shown in Fig. 5(a). The major contribution arises 
due to the hybridized states of Bi p and Sel p states 
situated on both sides of the first quintuple layer [^ . 
Interestingly, the contribution from the interface of the 
first and second quintuple layers is most significant. The 
Dirac states consisting of Sel (px^Py) and Bi pz states 
provide the largest contribution indicating their surface 
nature that is linked to the bulk via hybridization. In 
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addition, there is significant contribution from Sel p^-Bi 
{PxiPy) hybridized states. It is evident that the Dirac 
states in this system are not purely surface states, they 
are contributed by both surface and interface states in¬ 
volving top quintuple layer. 

Ironically, the details of the Bi terminated case de¬ 
picted in Fig. 5(b) exhibit an unusual scenario of Dirac 
states arising essentially from the interface states of first 
and second quintuple layers. This is different from the 
predicted surface behavior and demonstrate the com¬ 
plex nature of real materials. Here, the Sel Pz states at 
the interface are strongly hybridized and form the Dirac 
states. In both the above cases, the Se layers at the inter¬ 
face of the successive quintuple layers appear to be quite 
strongly coupled. This explains why cleaving of these 
materials require strong force, thereby, employment of 
top-post removal method. 

The presence of a layer of oxygen on Bi-terminated 
surface leads to significant change in the characteristics 
of the Dirac states. The Dirac states exhibit dominant 
surface character as shown in Fig. 5(c) along with signif¬ 
icant reduction of the interface contributions. The large 
electronegativity of oxygen compared to all other con¬ 
stituent elements leads to an effective pulling of the elec¬ 
tron cloud toward the surface. These results unambigu¬ 
ously demonstrate a way to engineer topological order 
on real systems, which can remain stable for long and 
can be used for various device application, realization of 
novel physics etc. 

In summary, the density functional calculations of 
Bi 2 Se 3 exhibit varied electronic structure on differently 
terminated surfaces. Although energy bands for every 
surface termination indicate cases of strong topological 
order, the signature of Dirac cone appears on Sel, Se2 
and Bi terminated surfaces as defined by SI, S3 and S4 
cases. We discover that the Dirac states are contributed 
by both sides of the top quintuple later with dominant 
contribution arising from the interface states rather than 
surface states of the pristine sample. Oxygenation leads 
to a shift of the Dirac point to higher binding energies 
for Se-terminated surface and makes the topological or¬ 
der weak. However, the Bi terminated surface is most 
robust and evolves to a better defined topological order. 
Robustness of topological surface states are the key fac¬ 
tor for the technological applications of these materials 
and realization of various exotic physics such as Major ana 
Fermions, magnetic monopole etc. These results demon¬ 
strate ways of engineering Dirac states on real materials. 


which lays the all important building block in this field. 
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